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Abstract 



O 
CT; 

5^ I We analyse exhaustively the structure of non-degenerate Cauchy 

horizons in Gowdy space-times, and we establish existence of a large 
class of non-polarized Gowdy space-times with such horizons. 

Added in proof: Our results here, together with deep new results 
of H. Ringstrom (talk at the Miami Waves conference, January 2004), 
■ establish strong cosmic censorship in (toroidal) Gowdy space-times. 



1 Introduction 



In 1981 Vince Moncrief pointed out the interest of studying the Gowdy 
metrics as a toy model in mathematical general relativity, and proved the 
fundamental global existence result for those metrics [15]. He developed 
approximate methods to study their dynamics, and discovered the leading 

'Partially supported by a Polish Research Committee grant KBN 2 P03B 073 24; email 
piotrOgargan. math. univ- tours .f r, URL www .phys .univ- tours . f rX^piotr 

^Supported by a grant from the Natural Sciences and Engineering Research Council of 
Canada; email lake@astro.queensu.ca . 



1 



order behavior of a large class of solutions of the associated equations [11]. 
Together with Beverly Berger he initiated the numerical investigation of 
those models [2]. In spite of considerable effort by many researchers, the 
global properties of those models are far from being understood. It is a 
pleasure to dedicate to him this contribution to the topic. 

In Gowdy space times the essential part of the Einstein equations reduces 
to a nonlinear wave-map-type system of equations [10] for a map x from 
{M,gap) to the hyperbolic plane {Jif,hab), where M = [T, 0) x with 
the flat metric g = —dt'^ + dO"^. The solutions are critical points of the 
Lagrangean 

^[A = \l tg^^'^habdax'^dfsx^dedt. (1.1) 

It is sometimes convenient to use coordinates P,Q G M on the hyperbolic 
plane, in which the hyperbolic metric hab takes the form 

h = dP^ + e^^dQ\ (1.2) 

Let Xf = ^,X0 = Z) denote the Lcvi-Civita connection of hah, and 
Dg = ^ := Dxg, Dt = := Dxf The Eulcr-Lagrange equations for (1.1) 
take the form 

DXt DXg _ Xt 

~Dt ~ ^^^^ 

or, in coordinates, 

+ Vl^ o xd^.x^d^'x" = 

where the F's are the Christoffel symbols of hah-, and U = — d^. 

As already mentioned, V. Moncrief established global existence of smooth 
solutions on (— oo, 0) of the Cauchy problem for (1.3) [15]. This implies [6, 15] 
that a maximal globally hyperbolic Gowdy space-time (^, g) with toroidal 
Cauchy surfaces can be covered by a global coordinate system (t,6,x^) G 
{—00,0) X X X , in which the metric takes the following form-*^ 

g = e-^/^\t\-^^^{-dt^ + d0^) + \t\e^{dx^f + 2\t\e^Qdx^dx^ 

+\t\{ePQ^ + e-P){dx'')\ (1.4) 

^ The form (1.4) has been claimed by Gowdy [10] under the hypothesis of two commut- 
ing Killing vectors, with spatial topology. This is not quite correct as there are a few 
further global constants involved [6] even in the current case of vanishing twist. However 
those constants can be eliminated after passing from the torus to its universal cover 

, or when working in local coordinates, and this suffices for the local considerations of 
the proofs given here. 
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where the function 7 solves the equations 

da = -t{\Xt\^ + \X0f) , dei = -2h{tXt,Xe) . (i.s) 

The main question of interest is the eurvature blow-up - or lack thereof - 
at the boundary t = of the associated space-time. An exhaustive analysis 
of this has been carried out in [9] for the so-called polarised case, where 
the image of the map x is contained in a geodesic in the hyperbolic space. 
There exist only two results in the literature which prove curvature blow up 
without the polarisation condition in a Cauchy problem context: The first 
is due to one of the current authors (PTC), who proves uniform curvature 
blow-up [7] under the condition that the solution at t = to satisfies^ 

supi^(|Xi|2 + |X,|2)(to,e) < ^ . (1.6) 

The second one^ is diic to Ringstrom [22], who assumes, at i = io, smallness 
of the derivatives of Q together with a bound 

0<7?<|io|J't(to,^)<l-??, (1.7) 

for some strictly positive constant 77. He further requires |to| to be sufficiently 
small. 

The purpose of this paper is to study Cauchy horizons in Gowdy space- 
times. Recall that existence of Cauchy horizons is precisely what one wants 
to avoid in order to maintain predictability of the Einstein equations. We 
wish therefore to describe properties of those Gowdy space-times which pos- 
sess Cauchy horizons, as a step towards proving that generic initial data for 
those space-times will not lead to formation of Cauchy horizons. Our re- 
sults here, together with those in [4,14,20-22,24], give strong indications 
that this will be the case, though no definitive statement is available so far. 

To put our results in proper context, we start by recalling some results 
from [7] concerning the geometric properties of the associated maximal glob- 
ally hyperbolic solution {^,^g) of the Einstein equations. The following 
gives a meaning to the statement that the set {t = 0} can be thought as a 
(perhaps singular) boundary of the space-time: 

^The threshold 6"^/^ in (1.6) has been recently raised to 1/2 in [4,24]. 

^There is no explicit statement about curvature blow-up in [22]; however, this fol- 
lows immediately from the results in [22] together with the calculations in the proof of 
Theorem 3.5.1 of [7]. The results in [22] have been strengthened and generalised in [24]. 
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Figure 1: The truncated domains of dependence Cl^{'^). 

Proposition 1.1 Every future directed inextendible causal curve {a,b) 3 
s r(s) in reaches the boundary t = in finite proper time. 

Further the limit 

lime(r(s)) , 
s— >6 

where 6 is the coordinate appearing in (1.4), exists. 

In order to continue, we need some terminology. For to < let the set 
Cf^{'tp) be defined as (compare Figure 1) 

CO (V') = {to<t<o, -\t\ <e-^< \t\} . (1.8) 

We shall say that limj^o {^^f = Oi if 

lim sup \f{t,e)-a\=Q . (1.9) 

Such limits look a little awkward at first sight; however, they arise naturally 
when considering the behavior of the geometry along causal curves with 
endpoints on the boundary i = 0. In any case, the existence of such limits 
can often be established in the problem at hand [4, 7] . 

The following relates properties of the Gowdy map x to curvature blow- 
up in space-time (the result follows immediately from the arguments in the 
proofs of Theorem 3.5.1 and Proposition 3.5.2 in [7]): 

Proposition 1.2 Let -ip e be such that 

hm^,, (,:)|t2z?,Xe| = \imco^(^^^\t''D0Xt\ = limc;o^(^) = . (1.10) 
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// there exists 1 7^ v{^p) G M such that 

limcOj(V>)l^^*l = '"W ' (1-11) 

then the curvature scalar C^^^^C"^'^'' blows up on every inextendihle causal 
curve in with end point in {0} x {tp} x x . 



Remark 1.3 Condition (1.11) can be weakened to 

either \imsup \tXt\ < 1 or liminf \tXt\ > 1 , 

where the Umits sup and inf are understood in a way analogous to (1.9). 

Remark 1.4 Suppose that there exists a sequence of points (tj, 9,i) € Cf|,(V') 

with ti ^ and \tXt\{ti,9i) — 1 > e for some e. If (1.10) holds along the 

same sequence, one then obtains curvature blow up along this sequence, 
which suffices to obtain inextendibility of the space-time at (0, ';/')• Prom 
this point of view the essential conditions are thus (1.10), while the existence 

of the velocity function v{iIj) defined by (1.11) is actually irrelevant in the 
following sense: v matters only at points at which it exists and equals one. 

Recall, now, that we are interested in conditions which would guarantee 
that Cauchy horizons do not occur. Consider, thus, a Gowdy space-time 
with a Cauchy horizon 7i. It is then easily seen (see Proposition 3.1 below) 
that there exists a point ip such that all curvature invariants remain bounded 
along the timelike curves t (t, ^/;, x^, x^). Proposition 1.2 suggests then 
that the asymptotic velocity v should equal one at ■0- This is, however, not 
entirely clear, because of the additional hypotheses (1.10) made. Now, one 
would like to have a sharp Cauchy horizon criterion without any unjustified 
restrictions. Here we prove the following related statement: 

Proposition 1.5 Consider a maximal globally hyperbolic Gowdy space-time 
which is smoothly extendible across a Cauchy horizon Ti,. Then there exists 
a point tp £ , coordinates {P,Q) on the hyperbolic space, and a number 
-foo(V') such that 

limc.o^(^^ {P{t, V) + In \t\) = Poo(V') , limco^wQ = • (1-12) 
Further there exists a sequence — ^ such that 

{\t\dtP){ti,^)^i. (1.13) 
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To continue, one would like to understand precisely the geometry of 
the set across which extensions are possible. Polarised solutions are known 
which have w = 1 on an interval [Oi^Or], and which are extendible across 
{0} X (9i,6r) X X (cf., e.g., [8]). This begs the question of existence 
of non-polarised solutions with t> = 1 on an interval [6i,9r]- In Section 2 we 
provide a simple construction of non-polarised Gowdy metrics, extendible 
across such a smooth Cauchy horizon; this is rather similar to a construction 
of Moncrief [18] (compare [3]). 

Based on those examples, one would naively expect that extendibility 
always requires the existence of such an interval [9i,6r]- Further, in those 
examples, one has much better control of the geometry than (1.10); one 
has AVTD^''^'' behavior near {0} x [9i,6r], as defined at the beginning of 
Section 3. Again one could hope that AVTDto''^^ behavior always holds 
across Cauchy horizons, but this has not been proved so far. The main 
result of our paper is the following theorem, which partially settles the issues 
raised so far: 

Theorem 1.6 Consider a smooth Gowdy space-time that is extendible across 
a Cauchy horizon H. Then: 

(i) TC is a Killing horizon; i.e., there exists a Killing vector field which is 
tangent to the null geodesies threading H. 

(ii) Ti is non- degenerate if and only if there exists an interval [9i, 9,.] C [a, b] 
such that for 9 G [9i,9r] the velocity function v{9) exists and is equal 
to 1 there, with 

\tPt - 1| + \tP0\ + \te^Qt\ + \tePQ6\ < Cf^" , (1-14) 

for some C, ap > 0, in a certain coordinate system (P, Q) on hyperbolic 
space. 

(iii) In fact, when Ti. is non-degenerate there exists a {P, Q) coordinate 
system on the hyperbolic space in which the solution is AVTD^^''^^ 
near {0} x [9i,9r], with Qoo = over [9i,9r], and with the conclusions 
of Theorem 3. 2 below holding on [9i , 9r\ ■ 

Remark 1.7 Assuming smooth Cauchy data, we note that there exists feo € 
N such that if the metric is extendible across a non-degenerate then 
it is smoothly extendible there. An explicit estimate for ko can be found by 
chasing differentiability in the constructions in the proof below to ensure that 
(3.5) holds. Smoothness of the extension follows then from Theorem 3.2. 
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Similar results can be established for Gowdy metrics with (sufficiently high) 
finite differentiability, whether in a classical or in a Sobolev sense. 

Remark 1.8 Extendibility above is meant in the class of Lorentzian man- 
ifolds; no field equations are assumed to be satisfied by the extension. We 

note that a necessary condition, in the non-degenerate case, for extendibility 
in the class of vacuum metrics is analyticity of the metric. 

Theorem 1.6 shows that degenerate horizons H., if any, would correspond 
to a single point in the (t, 6) coordinates. For analytic metrics non-existence 
of degenerate horizons Ji in the Gowdy class of metrics follows from the 
arguments of the proof of Theorem 2 of [19], but the information provided 
by those considerations does not seem to suffice to exclude the possibility 
that K = for smooth but not necessarily analytic metrics. While we find 
it unlikely that degenerate horizons could exist here, we can exclude them 
only under some supplementary hypotheses on the behavior of the Gowdy 
map x: 

Proposition 1.9 Let ip be such that v{ip) = 1 and suppose that either 
a) there exists ap > such that we have 

\td0P\ + \te^dtQ\ < C\t\''^ (1.15) 

on (V'), or 
h) Equation (1-10) holds, or 
c) the solution is AVTof'^^ on Q° (V-)- 
Then: 

(i) If there exists a sequence ipi ^ ip such that v{ipi) 7^ 1, then there exists 
no extension of ^ across {0} x '(p x x S^. 

(ii) Moreover, if 

_(ii„,,))%e-«<«(^Q„(«)%0, 

then the curvature scalar Ca/s^ysC'^^'^^ blows up on every inextendible 
causal curve in (^, '^g) with accumulation point in {0} x {ip} xS^ xS^. 
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We note that the power-law condition (1.15) is justified for initial data 
near those corresponding to a fiat Kasner solution {Pq, Qq) = (— In \t\,0) by 
the results in [4]. 

This paper is organised as follows: In Section 2 we present our construc- 
tion of solutions with Cauchy horizons. In Section 3 we define AVTD^^'*^^ 
Gowdy maps, and we give proofs of the results described above. In Section 4 
we investigate further the blow-up of the Kretschmann scalar Caf^^sC'^^^^ ■ 
In Appendix A the Riemann tensor of the Gowdy metrics is given. 

Added in proof: It has been shown recently by H. Ringstrom [24] that 
(1.10) holds at all points 9 e for all solutions of Gowdy equations. The 
arguments of H. Ringstrom [24] also show that the collection of Gowdy initial 
data for which the set {9 : v{9) = ±1} has no interior is of second category. 
Consider any solution in that last class, and suppose it contains a Cauchy 
horizon. Proposition 1.9 shows that the horizon is non-degenerate, which is 
impossible by Theorem 1.6. It follows that generic (in the second category 
sense) Gowdy metrics on have no Cauchy horizons, which establishes 
strong cosmic censorship within this class of metrics. 

2 Two families of solutions with v = and with 

V = 1 

There arc by now at least two systematic ways [20, 22] of constructing 
reasonably general families of solutions with controlled asymptotic behavior 
satisfying (compare (1.7)) 

0<ri<\t\Pt{t,9)<l-r] 9 e S\ to<t <0 . (2.1) 

The proofs given there only cover situations in which^ \t\Pt is bounded away 
from zero. It turns out that there is a trivial way of constructing families of 
solutions for which tX^ approaches zero asymptotically, as follows: Let 

be any rotationally-symmetric solution of the standard wave-map equation 
from the three-dimensional Minkowski space-time (R^'^,77) into the two- 
dimensional hyperbolic space {J^,h). Set 

.r{l.O)=y{x° = 9,r = -t) . (2.2) 

^However, it seems that certain families of solutions with vanishing v can also be 
constructed using Fuchsian techniques (A. Rendall private communication). 
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Note that the replacement \t\ ^ r, 6 ^ brings the action (1.1) to that 

for a rotationally invariant wave-map as above, so it should be clear that x 
satisfies the Gowdy evolution equation (1.3). Now, the solution x will not be 
periodic in 9 in general (and it seems that there are no non-trivial smooth 
wave maps y at all which are periodic in x^). However, it is straightforward 
to construct periodic solutions for which (2.2) will hold on an open neigh- 
borhood of some interval {0} x [b, 27r - 6] C {0} x [0, 2tt] {0} x S^, with 
< 6 < TT, as follows: let x be as in (2.2), and for 6 G [6/2, 27r - 6/2] let 

f{e) = x{-b/2, 9) , g{9) = dtx{-h/2, 9) . 

Extend / and g to 27r-periodic functions if^g) in any way. Let x?^ be 

the solution of the Gowdy evolution equation (1.3) with initial data {f,g) 
at t = —6/2. Uniqueness in domains of dependence of solutions of (1.3) 
shows that x^ . will coincide with x in the domain of dependence of the set 
{t = -6/2 , 9 G [6/2, 27r - 6/2]}: 

x^- - = X on {-6/2 <t<0, h + t < 9 <2tt -h-t} . 

Clearly one can also construct solutions which will satisfy (2.2) on the union 
of any finite number of disjoint domains of dependence as above. 

The wave maps y can be obtained by prescribing arbitrary rotationally 
invariant, say smooth, Cauchy data ?/(0, r) and dty{0,r), and the solutions 
are always global [5]. In this way one obtains a family of solutions x of the 
Gowdy equations paranicterised, locally, by four free smooth functions (this, 
by the way, shows that the "function counting method" might be rather 
misleading, as this family of solutions does certainly not form an open set 
in the set of all solutions). The y's are globally smooth both in the t and 
9 variables, which implies that the x's display the usual "asymptotically 
velocity dominated" behavior with 

limiXt = limXt = , (2.3) 

the convergence being uniform in 9 for 9 & [a,b]. In particular one has 
curvature blow-up for causal curves ending on {0} x [a, b] x x by 
Proposition 1.2. 

Equation (2.3) shows that x has zero velocity in the relevant range of 0's. 
It follows from the results proved in [4] that all solutions which have v = 
on an interval I are obtained from the procedure above on an neighborhood 
of {0} x I. 
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The solutions just described can be used to construct solutions with 
t; = 1 on open intervals. (Here and elsewhere v is always the quantity 
defined as in (1.11), whenever it exists, with Xf being associated with the 
{P,Q) representation of the solution as in (1.4).) Indeed, let x = {P,Q) be 
one of the zero- velocity solutions defined by (2.2), and define a new solution 
X by performing the "Gowdy-to-Ernst" transformation [21]: 

P :=-P- In \t\ , e^dtQ := -e^deQ , e^deQ := -e^dtQ . (2.4) 

The new map satisfies again the Gowdy equation (1.3). It immediately 
follows that we have 



\tXt\ - 1 +\tX0\< c\t\ 



so that we obtain power-law blow-up, together with the new asymptotic 
velocity equal to one on any interval on which the old velocity v = 0. All 
the resulting space-times have a smooth Cauchy horizon across any interval 
on which v = ±1, which can be checked by standard calculations using (3.6) 
(compare [8, 16 18]). 

As already mentioned in the introduction, a similar techniqTic, mapping 
singular solutions to ones with Cauchy horizons and vice-versa, has been 
used by V. Moncrief in [18], compare [3]. While V. Moncrief used this 
method to produce singular solutions out of ones with Cauchy horizons, our 
approach is exactly the reverse one. 



3 Proofs 

We start with some terminology. We will be mainly interested in solutions 
of (1.3) with the following behavior (for justification, see [4,9,14,20] and 
(4.9)-(4.8) below): 

p{t,e) = -vi{e)hi\t\+Poo{0) + o{i) , 

' |t|2^i(^)fvQW + o(i)) , o<vi{e)^N 
Q{t,e) = Qoo{e) + \ |tr(')(Qin(e)in|t| + ^Q(^) + o(i)) , o<v,{e)eN 

Qln(9)ln|i|+o(l) , „,(0)e_N*; 
. o(l) , -W^viie) < 



The function 



V := 



\vi\ 



(3.3) 
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will be called the velocity function, while Qoo will be called the position 
function. We shall say that a sohition is in the AVTO^^'*^^ class if (3.1)- 
(3.2) hold with functions vi, P^o, Qoo and iJjq which arc of Ck differentiability 
class (on closed intervals the derivatives are understood as one-sided ones 
at the end points). For A; > we will assume that the behavior (3.1)-(3.2) 
is preserved under differentiation in the following way: 

yo<i + j<k di{tdty(^p{t,e) + vi{e)in\t\~p^{e)j =o{i) , (3.4) 

similarly for Q. 

We continue with the 

Proof of Proposition 1.1: The fact that inextendible causal curves meet 
the set t = in finite time is proved in Proposition 3.5.1 in [7]. Next, since 
f is a time function on we can parameterize F by t: 

F(i) = (i,0(i),x«(t)). 

Timelikeness of F together with the form (1.4) of the metric gives 

d9 

which implies that 9 is uniformly Lipschitz along F, and the existence of the 
limit := limt^o 9[t) immediately follows. 

□ 



We are ready now to pass to the proof of the following result, mentioned 
in the Introduction: 

Proposition 3.1 Consider a Gowdy space-time with a Cauchy horizon H. 

Then there exists ip lE such that all curvature invariants remain bounded 
along the timelike curves t {t,^p,x^,x'^), —l<t<0. 

Proof: Let p G TC and let F be any timelike curve starting from p and 
entering the globally hyperbolic region We set F := Tf].^; by Proposi- 
tion 1.1 there exists & such that 9 approaches ^p along F. By hypothesis 
the extended metric is smooth around p, hence there exists a neighborhood 
^ of p on which all curvature invariants are bounded. In particular all cur- 
vature invariants are bounded on H Since the metric is U{1) x U{1) 
invariant on all the curvature invariants will also be bounded on the 
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set obtained by moving ^ PI ^ by isometries. This last set contains all the 
timelike curves as in the statement of the Proposition, with — e < t < 0, and 
the result easily follows. □ 



One of the ingredients of the proof of Theorem 1.6 is the following result, 
proved in [4]: 

Theorem 3.2 Suppose that there exist constants C, > such that 

\tdeP\ + \te^dtQ\ < C|t|"^ (3.5) 

and consider any point -0 £ [0',b] such that vi{ip) = 1. Then the functions 
{P,Q) := (P + In (5)|co can be extended to an AVTD^'^^ map from 

to Jf2- If limj^o = for all j G N, then for alli,k eN we have 

hm^o ^^)d^'+'dl{P + In \t\) = hm^o ^^^df+^d^eQ = ^^^C? W^e+'Q = « • 

Further, if vi = 1 on an interval [9i,6r], then the restriction {P,Q) := 
[P + In g^^y) can be extended to a smooth map from to J^, 
with (3.6) holding for all i/j G [6i,0r]- 

Remark 3.3 The vanishing of the last term in (3.6) for all /c > is some- 
what surprising. As already pointed out in the introduction, the power-law 
condition (3.5) is justified for initial data near those corresponding to a flat 
Kasner solution (Pq, Qq) = (— In 0) by the results in [4]. 

We are ready now to pass to the 

Proof of Theorem 1.6: Recall that H is a Kilhng horizon if there exists 
a Killing vector field which is tangent to the generators of TC. The fact that 
7^ is a Killing horizon follows from the proof of Proposition 1 of [1] (see 
in particular Lemma 1.1 and Lemma 1.2 there); further, any other Killing 
vector is spacclike on TC. Thus, if a Gowdy space-time is extendible across a 
connected Cauchy horizon, then one of the Killing vectors, say X, is tangent 
to the generators of the event horizon. Since our claims are local, without 
loss of generality we may assume that H is connected. Locally near a point 
in Ti we can construct a null Gauss coordinate system (T, W, x'^), A = 1,2, 
as in [19], leading to the following local form of the metric 

g = 2dT dW + Tcf) dW^ + 2T/3Adx^dW + HABdx^dx^ , (3.7) 
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with some smooth functions ^, Pa, and hab which have obvious tensorial 

properties with respect to x"^-coordinate-transformations. More precisely, 
let {W,x^) be local coordinates on Ti such that dw = X, and such that 82 
is another linearly independent Killing vector. Let k be the null vector field 
defined along H, pointing towards the globally hyperbolic region, such that 

g{dw,k) = l, g{dA,k) = 

on H. The local coordinates {W, x"^) are extended from 7^ to a neighbor- 
hood thereof by requiring {W, x^) to be constant along the geodesies issued 
from Ti. with initial velocity k. Letting T be the affine parameter along 
those geodesies, with T = on H, one obtains (3.7). Now, isometries map 
geodesies to geodesies and preserve affine parameterisations, which easily 
implies 

dw9fiu = = 829^,^ , 

throughout the domain of definition of the coordinates. Equation (2.9) 
of [19] for R^b shows that on any connected component of H there exists a 
constant k>0 such that 

(f)\n = K. 

Whatever the range of the W and coordinate in the original extension, 
we can without loss of generality assume that the functions above are defined 
for W, in M - or in - since those functions are independent of W and 

anyway. There might be difficulties if we are trying to build a manifold 
by patching together the resulting local coordinates, but we do not need 
to patch things back together, so this is irrelevant for the local calculations 
that follow. 

For r > we replace the coordinates T and W by new coordinates 
(t,x^), t > 0, defined as 

W = x^ + alni , T = /3P , (3.8) 

with constants a G M, /3 > 0, leading to 

g = a(3{'i + Ka + 0{P))dP + 2p{2 + aK + 0{P))tdidx^ 
+2i/3Adx^{tdx^ + adi) + PP{k + 0{t^)){dx^f 
+HABdx^dx^ . (3.9) 

We can choose a, /3 so that 

a/3(4 + Ka) < , 
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and we assume that some such choice has been made. 
Let Xi = dw = d^z , X2 = 8^2 , and define 

Kb = 9{Xa,Xb), a, 6 = 1,2. 

Prom (3.9) we have 

Xu = + Oii^) , X12 = PlP + 0{i^) , A22 = /X22, (3.10) 

Note that X2 is spacclikc on TC. It follows that 



(Vdet a) 



0{t), 
I 0, 



df, = df, K / 0, 

df, = di, K = 0, 

df, = d^i, Kj^ 0, 

da 



(3.11) 



d^i, K = 0, 
dj.i, i = 2,3, 



and 



t :-- 



Vdet A 



(3.12) 



0{P), K = 0. 

Now, away from the set t = the space-time metric g can be written^ in 
the form (1.4). The functions t and 6 are defined uniquely up to a single 
multiplicative constant, cf., e.g., [6]; the normalisation (3.12) gets rid of that 
freedom. By a rotation of the Killing vectors d^a we can always achieve 

d^i = X\ = , d^2 = X2 = dj.2 . 

Recall that a Killing horizon is said to be degenerate if (giX, X)) = 0. 
Since ^22 does not vanish at t = 0, (3.12) shows that TC is degenerate if and 
only if K vanishes. From now on we assume that k 7^ 0. Inspection of (3.9) 
shows that a convenient choice of a and /3 is 



2 + an = , a/3(4 + na) 



-1 



so that (3 = k/4. Ordering the entries as {i,x^,x^,x'^), (3.9) takes the 
following matrix form 

/ -i + o(P) 0(^3) oii) o{i) \ 

o{i^) pKp + oii^) oiP) o{P) 
o{t) o{F) Mil m 

V 0{i) 0(f2) /X22 / 
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This gives det g = —/3kP det /j, + 0{i^) and 





-1 + 0{P) 


0{i) 


0{t) 


o(i) 


\ 




o{i) 


{^kP)-^ + 0(1) 


0(1) 


0(1) 






o{i) 


0(1) 








V 


0(i) 


0(1) 









(3.13) 

where /x^^ is the matrix inverse to fiAB- Prom (3.13) and (3.11) we find for 
t < 



^ X ft /dt\ „ f.'idt dt 



dt 



. (-i.o(.))(|)%o(,|^.(.n,o(.))(- 

= -/?K;U22 + O(i^) • (3.14) 
Comparing (1.4) and (3.13) we also obtain 

/ \ f)ft Bf) 

= g{Vt, V9) = - [^fp^2 + OiP)) + Oii) ^ , 



and 



-givt,vt) = give,ve) 

di) didx^ 



de 



= -(i + 0(£2))(^|y + o(f) 



de ae 



89 Y 
dx^ J 



It follows that ddjdx^ is uniformly bounded from above and away from zero 
for t > 0, with did = Oit). This implies that Q extends by continuity to 
a Lipschitz function on and also implies the existence of the claimed 
interval of ^'s. 

Comparing (3.10) with (1.4) we find 

\t\e^ = PkP + 0{i^) , \t\e^Q = Pif + O(t^) . 

For K / this gives 



l^i + 0{P) = ^ + Oi\tf), Q=^ + 0{P), 
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and what has been shown so far about t and 9 further gives 

tPt = l + 0{t^) , tPe = 0{\t\) , te^'Qt = 0{\tf) , te^'Qe = 0{\t\^) . 

We have thus obtained a representation of the solution for which 

vi = -l, v = l, ItXel = 0{\tf) , (3.15) 

and the sufficiency part of point (ii) is estabUshed. At this stage one can 
directly derive a full asymptotic expansion of the solution using the Gowdy 
equations (1.3); this will lead to AVTD^''^^ behavior with a perhaps non- 
constant function Qqo- An alternative way consists in swapping the order 
of the Killing vectors; (3.15) will then still hold except for a change of sign 
of Vl, so that the hypotheses of Theorem 3.2 are satisfied, and point (iii) 
follows. 

Consider, finally, a solution with w = 1 at a point ip, or on an interval 
[01, Or], with (1.14) holding there. By Theorem 3.2 the solution is AVTB^'"^^ 
on Cj'Jj('i/'), or on [to,0) x [6i,0r], and by integration of (1.5) one easily finds 

7 = -ln|i| + 0(l) . (3.16) 

As already pointed out, such solutions are smoothly extendible by the cal- 
culations in [8]; alternatively, one can run backwards the calculations done 
above. In any case (3.14) holds. Prom (1.4) we have 

5(Vi,Vt) = -eT/2|i|V2, 

and comparing (3.16) with (3.14) the non-degeneracy of the horizon follows. 

□ 



The reader will have noticed that the last part of the proof of Theo- 
rem 1.6 shows the following: 

Proposition 3.4 AVTd'^^''^^ space-times do not contain degenerate hori- 
zons. 

We continue with the 

Proof of Proposition 1.5: The result follows immediately from the 
calculations in the proof of Theorem 1.6, with the following modifications: 

we choose X2 to be the Killing vector which is tangent to the generators 
of H, and we let Xi be any other Killing vector. If ^ is a point as in 
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Proposition 3.1, then g{Xi,Xi) tends to a strictly positive limit along the 
curve r defined in the proof of Proposition 3.1. Further this limit is the 
same for any causal curve which accumulates at the point p of the proof of 
Proposition 3.1. It follows that along any of the curves t {t^ip^x"-) the 
limit \\m.i^Qg{Xi,Xi) exists, and in fact one has that the limit 

exists. Using (1.4), this can be translated into the statement that the func- 
tion te^^^''^^ has a finite limit as t goes to zero, which justifies the first 
equation in (1.12). Since X2 is normal to H the function g{Xi,X2) vanishes 
on H, which implies the second equation in (1.12). Finally if \t\dtP avoids 
the value one, then either < 1 — e or \t\dtP > 1 + e for some e > for 

t large enough, which is incompatible with the first equation in (1.12). □ 

We close this section with the 

Proof of Proposition 1.9: Without loss of generality we can choose the 
coordinates {P,Q) on hyperbolic space so that wi(V') = 1. It follows from 
the results in [4] that under any of the conditions of Proposition 1.9 the 
conclusions of Theorem 3.2 hold; thus x is AVTD^''^^ in Cli^). Under 
the hypotheses of point (i), suppose that the solution is extendible. Since 
the velocity function is not equal to one in an open interval containing ip, 
Theorem 1.6 shows that the horizon must be degenerate. This contradicts 
Proposition 3.4, and establishes point (i). 

Point (ii) can then be established by inspection of the curvature tensor, 
which we give for completeness in Appendix A. The relevant calculations 
have been done using Grtensor [13]. The interested reader will find MAPLE 
worksheets and input files on URL http://grtensor.phy.queensu.ca/ 
gowdy. It follows from the the AVTB^'*^^ character of the solution on 
^to (V') that there exist bounded functions 700 and Z such that 

7 = -ln|i|+7oo(^) + o(l) , 

P{t, e) = -v{e) In \t\ + P^{9) + tZ{t, 9) . 

The field equations further show that \mit^QtdtZ{t,^) = (in fact Z = 
0(|i| In |f|); see the next section for more detailed expansions), and one finds 

i- i> + ^'''-^^^ (wQoo wY) e^^*''^) 

Cap.ysC'^f'^'^A^ ^ , (3.17) 
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as t — 0, provided that the coefRcient in the biggest parenthesis in the nu- 
merator does not vanish. As e'*'^*'^) ~ the curvature scalar Cq^^^C"^''''' 
diverges then like (As already pointed out in the introduction, in 

the next section we will study in detail the remaining possible behaviors of 



4 The blow-up structure of the Kretschmann scalar 

Throughout this section we assume that V is such that ui(V') = 1- To 
analyse the behavior of CaP'ysC"'^'^^ in the case 

-(^„W))%e-»(*)(^O.w)^0, (4.1) 

the first question to answer is how far to push an expansion of P and Q 
to isolate the potentially unbounded terms in the curvature. A tedious 
but straightforward inspection of the curvature tensor, as given in the ap- 
pendix, shows that terms of the form f{9)t'^ In^ \t\ in Q, with j ^ 0, might 
lead to a logarithmic blow up of the tetrad components R (i) (3) (1) (4) and 
(2) (3) (2) (4) (here the tetrad given at the beginning of Appendix A is 
used), and that any higher power of t will lead to a vanishing contribution 
to the tetrad components of the Riemann tensor. This means that one needs 
to have the exact form of all the coefficients in an asymptotic expansion of Q 
up to order 0{t^). Similarly one finds that one needs to have the exact form 
of all the coefficients in an asymptotic expansion of P up to order 0{t^). 
Those expansion coefficients can be found by collecting all terms with the 
same powers of t and of In \t\, and setting the result to zero, in the Gowdy 
equations, 

d^P - dlP = + e^' {{dtQ? - {deQf) , 

a|Q -djQ = -^-2 {dtPdtQ - dePdeQ) . (4.2) 

Recall that we are dealing with AVTD^''^^ solutions. It is easily seen from 
(3.1)-(3.2) with fi(V') = 1 and from (4.2) that (compare (3.6)) 

(5eQoo)(V') = 0. (4.3) 

Assuming this together with (4.1), we note the following simple observations 
concerning the behavior of the solution at 9 = ip: 
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(i) The term does not give any undifferentiated contribution to P. 

(ii) A term f{6)\t\'^^^^ W \t\ in Q witii a larger than or equal to, say 7/4, 
generates in P, for v{0) close to 1, terms of the form 

^^2v{e)+a{e)-2 f^j^Q^ ^ Jq^p^ powers of In \t\ 

+ higher powers of \t\ (multiplied perhaps by higher powers of In \t\) . 

(4.4) 

(iii) A term f{d)\tf^^^ W \t\ in P with /? > generates in Q, for v{d) close 
to 1, terms of the form /(0)|tp"(^)+'^(^) In-'^^ as well as terms with 
the same power of \t\ but lower powers of ln|t|, or terms with higher 
powers of \t\. 

Consider the linear counterpart of (4.2), 

d!f-dlf = -^. (4.5) 

As shown in [12] (compare [9, Equation (4a)]), for every smooth solution 
of (4.5) on (— oc, 0) x there exist functions f\n{t'^,0),f{t^,9), which are 
smooth up to boundary on the set (t^, 9) G [0, oo) x S^, such that 

/(^,^) = /ln(^^e)ln|^| + /(^^e). 

Further for any /in(0, 6) and /(0, 6) there exists a solution as above, and we 
have the asymptotic expansion 

/ = /i„(0,())ln|fH-/(0,9) 

Under the current hypotheses the linearisation of the P equation differs from 
(4.5) by terms decaying sufficiently fast so that the leading order behavior 
of P is correctly reflected by the above. That is not the case anymore for 
Q, because of the 1/t behavior of the d^P term, so that the leading terms 
in the Q equation linearised with respect to Q are 

dlf-d'ef = +^. (4.6) 

The associated indicial exponents are zero and two, from which it is not too 
difficult to prove the following behavior of solutions of (4.6) 

f{t,e) = Mt\e) + Mt^9)tHn\t\ , 
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with freely prescribable functions fo{0,0) and dgfo{0,6), together with an 
associated asymptotic expansion 

/ = fo{0,0) + -^^—^ 1 ln\t\ + -^^—^ {0,e)t +0{t ln\t\) . 

Since the full solution (P, Q) is already known to belong to the AVTD^''^^ 
class, using what has been said one can proceed as follows: one starts with 
the leading order behavior of P, 

P = -vi{0) Int + Poo W + OinCltl') , (4.7) 
where / = Oin(|i|*) denotes a function which satisfies 

iwi<^MiiriMiii)r^^ 

for some constants Cij,Nij. Inserting (4.7) into the equation for Q in (4.2) 
one finds that aX 9 = ip we have the expansions (recall (4.1), (4.3)) 

+2 

Q{t,e) = Qoo{9) + jdlQoo{9)ln\t\+ijQ{e)t' + t''W{t,9). (4.8) 
Inserting (4.8) in the first equation in (4.2) one then obtains 

P{t,9) = -v{9)ln\t\+P^{9)-jd^v{9)ln\t\ 

^(in^ \t\ - 2 In It I + -y 



+ 

+e2^- W {^q{9) + dlQ^{9)) ^^^^ ( In \t\ - l) 

I \^ dlP^{e) + d>{e) ^ ^2P.(,) ^ e 

+t^Z{t,e). (4.9) 

All the functions appearing above which depend only upon 9 are smooth, 
with f = Z 01 W satisfying estimates of the form 

{tdtr{doyf = Oi\\n\t\\^+'') 

for some A^. One can now insert those expansions in the Riemann tensor 
and obtain its behavior for |t| small. In the case —^j^ 7^ but = 
g2Poo(V')(^g^(^^))2^ a Grtensor calculation gives, again at 9 = tp, 

C^p.sC"'"'' ~ 3t (ln(|t|))4e^(*'^) I^^^Q^ {9)^' (e^-(^))' (4.10) 
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so that with e^^*'^) ~ l/\t\ then Cap-ysC"^^^ now diverges hke (ln(|t|))^. 
Next, if u(V') = 1 and = ^Qoo(^) = we obtain 

Cap-.sC''^''' ~|t|(ln(|t|))2e^(*'^)i?(^) (4.11) 
as i — 0, provided that does not vanish. Here 

(4.12) 

Since e<^^^^^ ~ l/|t|, we find that C^p^sC"^^^ diverges hke (ln(|t|))2. 
One can hkewise check what happens if -R(V') = 0: 

Co.p.sC'^^'' ~ -3(^^(^))(2(^Poo(^))' + 2^Poo(V) +4e^-W^VQ(V') 

+8e^-WVQW^i'oo(V') + ^vm\t\H\t\y^^'''^^ , 

and one obtains ln|t| behavior of the Kretschmann scalar unless the coeffi- 
cient above vanishes; in that last case the Kretschmann scalar is bounded. 



A Appendix. The curvature tensor of Gowdy met- 
rics 

We use the following tetrad 

ei»= [e(V4 7(t,e))t(i/4)^ 0, 0, (| 
e2^= [0, e(V4 7(M))i(i/4)^ q, O] 



e5" = 



0, 0, -Q(t, 6) 



oP{t,e) 



0, 0, 







Writing a partial derivative as a subscript, a GRTensor calculation 
with Maple gives the following components of the Riemann tensor in this 
frame: 

Rii) (2) (1) ^2) = -/^/''){-{en't^Qt'-2tHen^Qt'Pt-2t^en^QtQt,t-t'Pt' 

+ {e^r Qt Qe Pe + 2t^ {e^f Qe,e) /t^''^^ 
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R{i) (2) (3) (4) = -^e^Vie(V27)(_p^g^ + Q^P^) 

Rii) (3) (1) i3) = -le(^/'^){-5tPt-At^Pt,t + tHen'Qt'Pt + t^Pt' + tHen'Qe'Pt 
+ 3t'Pt Pe^ + 5 {eP f Qt^ + Pt' - ^ {e^f Q,^ - P,H' - 1 + 2t^ {e^f Q, Qe Pe) 

(3) (1) (4) = Ve(V27)(_8tP,Q, + Q,3^2(gP)2 + g^^2p^2^3Q^^2(gP)2g^2 

+ Qtt^P0^ -5Qt-4tQt,t + 2t^QePePt) M 

R(i) (3) (2) ^3) = le(y'^\4tPt,e-St^Pt''Pe-2t-'Pt{eP)^QtQe-^t{ePfQtQe 

- Pe e {en' Qt^ - Pe ^ {e^f Qe' - Pe' t^ + 3Pe)/Vi 

R{i) (3) (2) (4) = -le''e('/^^\6tQtPe-2Qtt''PePt-3Qt't^en'Qe + 4.tQt,e 

- e Qe Pt' - Qe' {e^f - f Qe Pe' + SQe + 2Qet Pt) j^t 

RW (4) (1) ^4) = le(y^'^\l-5tPt-At'Pt,t + 3ien't^Qt'-t'Pt' + tHen'Q0' + Pe't^ 

+ t' (e^)2 + ^3 p^3 + ^3 (gP)2 Q^2 + 3 ^3 p2 + 3 ^3 (gP)2 g^ p^) ^ 

t(3/2) 

^(1) (4) (2) (3) = -^e^e(l/27)(2tQ,P,-2Qtt2p^p^_3g^2^2(gP)2g^ + 4^g^^^ 

- f Qe Pt' - f Qe' {e^' - 1' Qe Pe' + SQe + QQet Pt) /Vi 

R{1) (4) (2) (4) =-^e(V27)(4tp^ ,_3t2p2p^_2i2p^^gP)2g^g^_4^(gP)2g^g^ 

- Pe e (e^' Qt' - Pe t' {ePf Qe' - Pe"" + 3 Pe) /Vt 

i?(2) (3) (2) (3) = -^e(V27)(t3(gP)2g^2p^+^3p^3+^3(,P)2g^2p^ + 3^3p^p^2 

- (e^' t' Qt' - t' Pt' + 5*2 (e^)2 g,2 + p^2 ^2 _ ^p^ + ^ _ 4^2 p^_^ 
+ 2t' {e^'QtQePe) /t^'"^ 

i?(2) (3) (2) (4) = ^e^e(V27)(g,3^2(,P)2 + g^i2p^2 + 3g^^2(,P)2g^2 + g^^2p^2_g^ 

-StPeQe + 2e QePePt-^tQe,e)/Vi 
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R (2) (4) (2) (4) = ^e(V2T)((eP)2 ^2 q2 + ^2 p2 + 3^2 (,P)2 g^2 _ p^2 ^2 + ^3 (,P)2 g^2 
+ Pt' + (e^)2 + 3^3 p^2 _ 1 _ ^p^ _ 4^2 p^_^ + (e^)2 Q, Q, p,) 

1 e(V2 7) (^2 p^2 _ 1 ^ (^P)2 ^2 q2 _ p2 ^2 _ ^2 ^^P^2 g^2^ 
R (3) (4) (3) (4) = - 4 ^(57^^ 
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